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Abstract 

Hard diffractive electroproduction of longitudinally polarized vector mesons 
is calculated within the leading In -r^ — approximation of QCD using the 

"-QCD 

leading order parton densities within the nucleon. Novel QCD features of 
the production of excited states and of the restoration of flavor symmetry 
are analyzed. At the onset of the asymptotic regime, our analysis finds an 
important role of quark Fermi motion within the diffractively produced vec- 
tor mesons, and we suggest to use this effect to measure the high momentum 
tail of the wave function of the vector mesons. We deduce a kinematical 
boundary for the region of applicability of the decomposition of the hard 
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amplitudes over powers of and/or a limit on the increase of the cross 
sections of hard processes at small x, and briefly analyze its consequences. 

We also estimate the nuclear attenuation of the diffractivc clcctroproduc- 
tion of vector mesons and compare with estimates of the shadowing of the 
longitudinal structure function. 
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I. INTRODUCTION 



Diffractive (coherent) production of hadron states in deep inelastic lepton-nucleon scat- 
tering is a new kind of hard process calculable in QCD. It provides the unique possibility 
to study the properties of vacuum exchange in QCD as well as the interplay between soft 
(nonperturbative) and hard (perturbative) contributions to the Pomeron. One of the main 
aimes of such research is to obtain a three dimensional image of a hadron as compared to 
the one dimensional images extracted from the usual deep inelastic scattering off hadrons 
and nuclei. The idea to use diffractive electroproduction of vector mesons to investigate 
color coherence effects has been suggested many years ago in Refs. 

It has been shown in Ref. 0] that the cross sections of hard diffractive processes off 
a proton target are proportional to the square of the gluon distribution in the proton. 
Recently, Brodsky et al. M derived within the leading as In -jS — In - approximation a 
QCD prediction for the diffractive production of vector mesons built of light flavors in 
deep inelastic scattering. Diffractive photo- and electroproduction of J/^ mesons has 
been calculated by M.Ryskin |^ within the BFKL approximation of pQCD and within 
the charmonium constituent quark model approximation for the light-cone wave function 
of the J/\E' meson. 

On the experimental side, first HERA results on p-meson production aX Q"^ > 8 GeV^ 
0] and on photoproduction of J/\E' mesons ||^,|^ have confirmed the fast increase of the 
cross section with energy predicted by pQCD. 

In this paper, we extend the QCD analysis of Ref. 0] in several directions. We explain 
that the formulae of Ref. P| are valid within the conventional leading In — approx- 
imation. So the cross section for hard diffractive electroproduction of vector mesons is 
expressed through the conventional leading In -r^ — order gluon distribution, and not by 
means of its asymptotics within leading order in In «^lni as in Ref. IHl. This cal- 

^^QCD ^ 

culation justifies using the LO QCD improved parton model densities to compare QCD 
predictions for the electroproduction of vector mesons with experimental data. 

We analyze the average transverse sizes of the qq components effective in cr-y^p, Kl, and 
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in a^^p^vp, b-fi^Vi and find them to decrease strongly with increasing Q^. At Q"^ ~ 10 
GeV^, those sizes are comparatively close, h^j^ ~ &7£^v ~ 0.3 fm, while at larger the 
rate of decrease of this average size is larger for 6o-^. Thus we argue that the current 

~ 10 GeV^ data suggest that the measured radius of the color distribution within a 
meson is b^^^v ~ 0.3 ^ 0.4 fm. This observation justifies the application of pQCD 
for the calculation of hard diffractive electroproduction of longitudinally polarized vector 
mesons. Really, the screening of the color charges insures the smallness of the effective 
interaction due to asymptotic freedom. In the case of cr^^p, the effective transverse size 
of the qq components is expected to be intermediate between those in a^^p and those in 
photoproduction because the contribution of soft QCD is enhanced for cr^^p. Also, the 
significant difference between foo-^ and b^^^^v indicates that substantial next-to-leading or- 
der corrections should be present in diffractive vector meson photo- and electroproduction 
which, in turn, should lead to a change of the scale effective in the gluon density. 

We show that the QCD prediction for the dependence of the p-meson production 
cross section is substantially weaker than that given by the naive dimensional estimate 
of cr^*p^v'p(x, Q^) oc Q'^ as a result of the increase of the parton distributions due to 

evolution We find also a further slowing down of the dependence due to the 
transverse quark motion in the vector mesons. Account of both effects resolves, to a large 
extent, seeming contradictions of predictions of Ref. with the dependence of the 
cross section of the ZEUS data . Thus, investigation of the dependence of the cross 
section of electroproduction of vector mesons may turn out to be an effective method to 
probe the three dimensional distribution of color in hadrons. 

We analyze the energy and dependence of the restoration of SU (3) flavor symmetry 
in the electroproduction of vector mesons built of different flavors. We find that the ratio 
of the yields of (p- and p-meson production should increase with increasing and/or 
1/x. We explain also that, in hard diffractive processes, SU{A) flavor symmetry should be 
violated in the opposite way as compared to the low energy regime - production of heavy 
flavors should be enhanced at sufficiently large Q^. To analyze the onset of the restoration 
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of SU{4:) flavor symmetry, we evaluate diffractive photo- and electroproduction of J/\E' and 
mesons within the charmonium model approximation to the light-cone wave function 
of the J/\E' and \E'' mesons, as suggested in Ref. 0, and find that account of the quark 
Fermi motion leads to a significant additional suppression of J/\E' photoproduction by 
a factor ~ 1/8 as compared to the asymptotic estimate, and to a slowing down of the 
onset of the asymptotic regime. Comparison with the data of Ref. ||TO|JTT| J^JT^ on J/\E' 
photoproduction shows that the pQCD predictions can be brought to agree with the data 
within the uncertainties in the gluon distribution of a nucleon. We explain also that the 
charmonium models differ from QCD for large Q^. We furthermore present estimates for 
the yields of excited vector meson states in diffractive processes, 'jI + p —>■ V + p. The 
interesting new effect, specific for QCD, is that the relative yield of excited states should 
increase with increasing ^ and/or Q^. 

Towards the end of the paper, we estimate the boundary of the kinematical region 
where, at very small x, the fast increase of the cross section with decreasing x should 
slow down and where existing methods of hard pQCD seem to be insufficient. We deduce 
the constraint on the kinematical range of applicability of the leading logarithm approx- 
imations and/or on the increase of the cross sections of hard processes with energy. We 
discuss practical consequences of this constraint and present model estimates for nuclear 
shadowing of the gluon distribution and of the diffractive vector meson electroproduction 
in interactions with heavy nuclei. These effects should be significant in kinematics which 
can be probed at colliding electron-nuclear beams at HERA and in the LHC heavy ion 
experiments. 

Diffractive photo- and electroproduction of vector mesons has been investigated re- 
cently in Refs. [p!3|-|T7[| within the constituent quark model based on the two-gluon- 



exchange approximation . Some of the predictions of that model - like the decrease of 
the effective size of the vector mesons with increasing or a significant production 

of p' mesons - resemble results obtained in this paper. Other predictions - as the dip 



in the p' meson production at small Q suggested in Ref. or the qualitative difference 
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between the behavior of radial and orbital (e.g. D-state) excitations - do not follow from 
the QCD formulae of Ref. |Q and the current paper. 

At the same time, we want to emphasize that at sufficiently large Q^, where a QCD 
calculation can be substantiated, the cross section should be expressed through the distri- 
bution of bare quarks in the vector meson and the parton distributions within the target ^ 
but not through the distribution of constituent quarks as in the models of Refs. |]5|,p!3Hl7|. 
There exists no direct relationship between the minimal Fock component of the light-cone 
hadron wave function, which enters in the QCD formulae for hard processes, and the wave 
functions of the constituent quark model. For example, the minimal Fock component of 
the quark-gluon wave function of vector mesons should decrease at large kf as a power of 

(cf. Ref. |]TP|) but not like a Gaussian as used in Refs. [P"^-[T7|. Just this difference is 
responsible for relevant effects found in the current work which are due to the transverse 
quark Fermi motion in the diffractively produced vector mesons. The qualitative differ- 
ence between QCD calculations and the two-gluon-exchange model of Ref. []TB[ has been 
discussed in Ref. [Q, and additional references to other works in that realm can also be 
found there. 

The organization of this work is as follows: In Sect. II, we give leading a^ln-r^ — 
predictions for the diffractive vector meson production in deep inelastic scattering, and, 
in Sect. Ill, we justify our pQCD approach listing the necessary approximations. In Sect. 
IV, we give estimates of the effects of transverse quark motion within the diffractively 
produced vector mesons and we discuss the NLO rescaling. In Sect. V, we compare 
our results with the few currently available data. In Sect. VI, we discuss the restoration 
of flavor symmetry and, in Sect. VII, the production of excited states. In Sect. VIII, we 
estimate a kinematical constraint on the region of applicability of the leading logarithm 
approximations, and, in Sect. IX, nuclear shadowing in the hard diffractive production of 
vector mesons as well as the longitudinal structure function of a nucleus are considered. 
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II. LEADING as In PREDICTIONS FOR THE DIFFRACTIVE VECTOR 

^QCD 

MESON PRODUCTION IN DIS 



It was demonstrated in Ref . @] , that a new kind of hard process, the coherent electro- 
production of vector mesons off a target T, 

7*+T^y + T, (1) 

is calculable in QCD for ^ My within the leading In -r^ — In - approximation. 

^^QCD ^ 

Here, V denotes any vector meson {p,u!,(j), J/"^) or its excited states. For comparison 
with previous works and appropriate references see Ref. [Q. 

The idea behind the calculation of hard diffractive processes is that, when the coherence 
length, Ic = 2mi^x ' exceeds the diameter of the target, 2rT, the virtual photon transforms 
into a hadron component well before reaching the target and the final vector meson V 
is formed well past the target. The hadronic configuration of the final state is then the 
result of a coherent superposition of all those hadronic fluctuations, |n), in the photon 
wave function whose masses satisfy the condition: 



« 1 • (2) 



3 

Here, tmin is the minimal momentum transferred to the target: 

Thus, as in the more familiar leading twist deep inelastic processes, the calculation 
should take into account all possible diffractively produced intermediate hadronic states 
satisfying Eq. (H). The use of completeness over those states allows us to express the result 
in terms of quark and gluon distributions, as in the case of other hard processes. The 
matrix element of electroproduction of a vector meson can thus be written as a convolution 
of the light-cone wave function of the photon, ip^^^n), the scattering amplitude of the 
hadron state, A{nT), and the wave function of the vector meson, ipy: 

A = A{nT) . (4) 
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In the case of a longitudinally polarized photon with sufficiently large Q^, the sum over 
intermediate states \n) is a qq pair. This can be demonstrated by transition into impact 
parameter space and by showing that the essential distances, b, between the quarks in 
the wave function of the photon are 6^ oc l/Q^. The situation is qualitatively different, 
however, for a transversely polarized photon due to the singular behavior of the vertex 
7t ~^ when one of the partons carries a small fraction of the photon momentum. In 
this case, soft and hard physics compete in a wide range of and x. 

Within the leading logarithmic approximation, i.e., as In -r? — In - ~ 1, when terms of 

^^QCD ^ 

the order as In -r^ — are neglected, the final result [H for the cross section for the produc- 

^^QCD 

tion of longitudinally polarized vector meson states, when the momentum transferred to 
the target t tends to zero, isQ 

_ 12n^ry^,+,-Mva'i{Q)r]l \ (l + ^) xGt{x, Q^) " 



dt 



(5) 



Here, Fy^e+e- is the decay width of the vector meson into an e^e~ pair. The parameter 
riv is defined as 

2jdzd^kt^viz,kt) ' ^ ' 

where $y is the light-cone wave function of the qq component of the vector meson. 

The factor \as{Q'^)xGT{x,Q'^)\'^ in Eq. (||) originates from the square of the cross 
section for the high-energy interaction of a small-size qq configuration with any target. 
It can be unambiguously calculated in QCD for low x processes by applying the QCD 
factorization theorem. In the approximation where only the leading as In — terms are 



accounted for, the result is 21,E2^ 



^In Ref. Q, a factor 4 was missed in the numerator of Eq. (^) |^0|. We are indepted to Z. Chan 
for pointing this out (A.Mueller private communication). 
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where b is the transverse distance between the quark q and the antiquark q, xGt{x,Q'^) 
is the gluon distribution in the target, and A ~ 10 according to our estimate in the next 
section. In this equation, the evolution and small x physics are properly taken into 
account through the gluon distribution. o"(6^) can be inferred from the calculations of Ref. 
||23| , which were performed in the leading a^lnx approximation of pQCD using certain 
input from hadronic quark models. Also, a quantity which can be related to cr(6^) has been 
calculated in Ref. within the leading a^lnx approximation of pQCD while making 



specific assumptions on the parton structure of the constituent quarks. 

The upper bound of the x-range of applicability of Eq. (^) is dictated by the condition 
that the formation time parameter should be much greater than the interaction time, 
which corresponds to the requirement that the longitudinal coherence length is larger 
than the target's diameter: 

1 

X < ^ 0.06 . 8 

Eqs. @ and can be easily generalized to the somewhat less restrictive leading 
as In ^^^a approximation. Derivation of the respective formulae is relatively simple for 
hard diffractive processes initiated by longitudinally polarized photons because longitudi- 
nally polarized photons produce quarks in a PLC (point-like configuration) only. It can 
be demonstrated by transition into impact parameter space that the essential distances, 
b, between the quarks in the wave function of a longitudinally polarized photon with 
virtuality are 6^ oc Using b as the small parameter, it is easy to determine 

which Feynman diagrams dominate for the hard diffractive process. It is important to 
note that the integral over the distance between the quarks produced by the longitudi- 
nally polarized photon does not lead to terms proportional to In 6^ in the evaluation of 
the box diagrams. So the whole contribution is concentrated at small b. As a result, the 
factorization/decoupling theorem of QCD can be applied to factorize the hard part of the 
amplitude. So, to calculate the amplitude in the leading In -rr— approximation, we 
need to accurately evaluate the hard blobe and then convolute it with the parton distri- 
butions. The hard blobe is described by the Feynman diagrams in Figs, la and lb, where 



^luons are attached to the quark box diagram. 




(a) (b) 

FIG. 1. Feynman diagrams relevant for the evaluation of the cross section of diffractive 
electroproduction of longitudinally polarized vector mesons, i.e., the jf^ + T^V + T process, 
in leading Os In — approximation. 



^QCD 



Direct calculation of the sum of Feynman diagrams for the 72 — > transition of Figs, 
la and lb where only two gluons couple to the hard blobe - although all interactions 
of these two gluons with the target are included - leads to seemingly the same formulae 
as in Ref. 0]. Really it is easy to demonstrate that within the leading a^lnfo^AQj^^ 
approximation, which is equivalent to the leading as In -j-^- — approximation, the dominant 
contribution is given by the diagrams with exchange of longitudinally polarized gluons. 
The diagrams corresponding to the exchange of transversely polarized gluons and those 
due to the quark sea in the nucleon are suppressed at small x. For a detailed discussion 
see Ref. . 

The major difference from 0] is that, within the approximation discussed above, the 
full gluon distribution in the target calculated within the leading as In -r^ — approximation 
should be used. In the formulae of Ref. |Q, however, xGriXjQ"^) is the gluon structure 
function calculated within the leading as In -jS — In - approximation. Remember, however, 
that corrections of powers are numerically small and negative Q]. Here, / is the gluon's 
transverse momentum within the "pomeron". 

In principle, the nondiagonal gluon distribution enters into those formulae. It will 
be explained in the next section that the usual parton distribution in the target T is a 
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reasonable approximation to this nondiagonal density. For a better understanding of the 
definition of the nondiagonal parton densities, these formulae should be compared to a 
similar one for the total photoabsorption cross section off a target T: 

cr^*T = / daGT{a,Q'^)a^*g{ai',Q'^) . (9) 

Jx 

Eq. (PI) is valid within the leading In -tS — approximation, and the Q"^ evolution is 
contained in the dependence of the parton distribution in the target T. 

The analysis of the Feynman diagrams for ImA^^j'^vT shows that, for diffractive 
vector meson production, the region of integration over a is restricted as: 

a = x + ^. (10) 

Here, 

^ ml + 

z{\-z) ^ ' 

is the invariant mass of the qq pair, and z and kt define the light-cone momentum carried 
by the quark/antiquark in the wave function of the 72- The second term in Eq. (|T^) can 
be neglected for the production of vector mesons built of light quarks, since effectively 
k\ <^Q'^ . This is because the wave function of the vector mesons decreases comparatively 
fastly with kl up to the onset of the pQCD regime. This explains why it is a reasonable 
approximation in JmA^^T^VT - but not in Eq. (|) - to pull the parton densities at a ~ x 
out of the integrand. 

The respective hard cross sections can be easily calculated through the miminal Fock 
components of the wave function of the vector meson. Thus, we deduce our asymptotic 
expression for the cross section of the electroproduction of vector mesons: 

2 

• (12) 



dt 



l2T:^Vy^,+,-Mval{Q)T{Q'')Til\{l + zf ^) xGt{x, Q' 



The factor T^Q"^) accounts for preasymptotic effects, i.e., T{Q^ — > oo) = 1, see the 
discussion in Sect. IV. 
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The major practical difference between the formulas deduced within the leading 
as In -rr — In - and the leading In — approximation (i.e., within the evolution equa- 

QCD ^ QCD 

tion) is that, in the leading ag In — approximation, xGt(x, Q^) should be interpreted 
as the conventional parton distribution. 

III. JUSTIFICATION OF THE PQCD APPROACH FOR DIFFRACTIVE 
VECTOR MESON ELECTROPRODUCTION IN DIS 

The formal derivation of the pQCD expressions for the cross sections of hard diffractive 
processes in terms of the parton distributions in the target consists of several steps. The 
first is to demonstrate that, at small x, it is legitimate to use completeness over hadronic 
states to express the amplitude in terms of hard interactions of quarks and gluons. We 
discussed this in the previous subsection. The second step is to prove that the process 
is really hard and that it is legitimate to apply the QCD factorization theorem, and 
the third step is to show that the nondiagonal parton distributions of the target can be 
approximated through the conventional parton densities. 

A. Applicability of pQCD and small transverse distances 

To better understand the issue of applicability of pQCD to the process (|I|), it is conve- 
nient to perform the Fourier transform of the amplitude into the impact parameter space 
of the qq pair, which leads to 

•Ay'^N-^VN oc J d^b dz i)l^*^{z,h) a{b'^) ipv{z,h) . (13) 

Here, the Sudakov variable z denotes the fraction of the photon momentum carried by 
one of the quarks, h is the transverse distance between the quark and antiquark within 
the photon, and o-(6^) = ^b'^as{b'^)xGN{x,P) is the color-dipole cross section of Eq. (0) 
off a nucleon target. The advantages of the description of high energy processes in terms 
of light-cone wave functions of bound states in the impact parameter representation has 
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been understood long ago (cf. H.Cheng and T.T.Wu and references therein). The 
qq component of the hght-cone wave function of a longitudinally polarized virtual photon 
reads 



(z, 6) = 2 Q z(l - z) Ko [Qh^z{l - z)j (14) 

in the limit of vanishing quark mass, and where Kq is the Hankel function of an imaginary 
argument. 

To estimate which values of h dominate in the integral, in Eq. (|1^) we approximate 

ipv{z,h) by 

tljv{z,h) (X z{l — z) jib Ki{jib) ^ (15) 

This form accounts for the common wisdom of the z dependence of the wave function of a 
longitudinally polarized vector meson, while the kt dependence is chosen^to be oc (pq^^;^- 
The parameter is related to the average transverse momentum of a quark/antiquark 



within the qq component of the meson's wave function through J {k^) = jJ,/V2. In our 



numerical analysis, we vary J {k^) between 300 and 600 MeV/c. Our guess is that, in 



reality, y {kf) should be of the order of 400 -r- 500 MeV/c, which is somewhat larger than 
the partons' average transverse momenta in hadrons. This is because, from the whole set 
of Fock components of the hadron wave function, the overlap with the wave function of 
the 72 emphasizes the component which has no soft gluon field. An analysis of the energy 
denominators relevant for the wave functions of vector mesons shows that, in order to 
avoid a spurious pole, the qq component's energy denominator may dominate only in the 
kinematics 

777-,, ~|~ kj- Q 

-f \ » . 16 

z(l - z) 



^In QCD, the expected behavior of oc ^ at large kt corresponds to even smaller interquark 
distances in the wave function of the longitudinally polarized photon. 
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This condition is restrictive only for the minimal Fock space component. When the number 
of constituents is large, the energy denominator is large more or less automatically. 

Now we need as{b'^)xGN{x,b'^). We estimate it based on the impact parameter repre- 
sentation for aL{x,Q'^), 



(^l{.x,Q ) = j dhdz agg-Ar(x,6) , (17) 

the nucleon's total longitudinal photoabsorption cross section, where 



2 

1^ 



TT 



12_ A 



In the derivation of this equation, to simplify the estimates, we again pulled out the parton 
distributions in the mean point x. 



We perform the analysis in two steps. First we use Eq. (18) and neglect the dependence 
of as(&^) xGn{x., b"^) on 6^ to estimate an average b{Q) for this process. The latter we define 
as the value of b up to which we have to integrate to saturate 50% of the integral in Eq. 
(0). We then identify Giy{x,b) ^ Giy{x,b{Q)), determine a new average b{Q), and then 
iterate this procedure to self-consistency. The final result of b{Q) for <7l{x,Q'^) obtained 
in this manner is plotted solid curve in Fig. 2. 

Thus, approximately A ~ 9.2 or 

9.2^ 



Gjqix^ 6^) ~ Gjq 



for X ~ 10^^. For fixed Q^, the transverse size b{Q) slowly decreases with decreasing x: 
between x = 10~^ and x = 10~^ it drops by about 10%. The above reasoning accounts for 
some of the next to leading order asf i In -r^ — ) effects, based on the intuitive idea that 
the interaction is determined mostly by the spacial region occupied by color. Note that a 
numerically similar result was obtained by Nikolaev and Zakharov |2^ in a consideration 
of Fj^ and OFt / dlxiQ^ (A ~ Ao- ^ 10). However, in the derivation of the quantity A(j in 
Ref. [^, the corresponding integrals were actually dominated by large 6's for which it 
is not clear whether the expression for the color-dipole cross section of Eq. (0) is at all 
applicable. Besides, no x-dependence of their A^, was discussed in Ref. [pH] . 
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FIG. 2. Average transverse size of the qq components effective in A-y'^N—fVN of Eq. ( |l3|) and 



of Eq. (|l^). In our numerical calculations, we used the GRV94 parametrization |27| of the 
nucleon's gluon density. 

Then we use this value of b{Q) to calculate the average transverse distances, b^^^p{Q), 



relevant for longitudinal p-meson production. The results of such a calculation for y (k^) = 
300 and 600 MeV/c are presented also in Fig. 2. One can see that b^^^p is quite small. 



however its numerical value is rather sensitive to the average value of y (kt). We observe 
also that b^^^p > b^j^. This difference may be rather small for ^ 10 GeV^, but it 
becomes significant for ~ 100 GeV^. 

Note that the values we find for b^^^p are substantialy smaller than the typical dis- 
tances between the valence quarks in a meson, 2r^/ ~ 1.2 fm. So, color is sufficiently 
screened within the wave function of the produced vector meson, and asymptotic free- 
dom can be used to justify the applicability of pQCD for the calculation of diffractive 
electroproduction of p-mesons at > 10 GeV^. 

Our results for the elementary color-dipole cross section of a small size qq pair with a 
nucleon target, (JggN^x, b) of Eq. (|TB|), are depicted in Fig. 3 for the GRV94 parametriza- 



15 



tion of the nucleoli's gluon density and for various values of x. As expected, the 
cross section rises approximately quadratically with h and it increases dramatically with 
decreasing x. Note, however, that only for the smallest x, the cross section <JqqN{x,b) 
approaches normal hadronic values 20 mb). However, for such x, it reaches those al- 
ready at quite small values of b, which, in turn, hints a limit of applicability of the leading 
logarithm approximations discussed here. This will be discussed further in Sect. VIII. 



I 1 



I I 

b 



30 



20 - 



10 - 







X = 10" 



GRV94 




"T"- ■'f' ''• ' ''' I I I I I I I I I I I I I I I I I I I I I i_ 



0.1 0.2 0.3 0.4 0.5 0.6 
b [fm] 

FIG. 3. Color-dipole cross section, cjqqi^{x, h) of Eq. (|l8|), as a function of the transverse size 



of the qq pair for various values of x and for the GRV94 parametrization |27| of the nucleon's 
gluon density. 



Results for an analogous evaluation of the average transverse distances relevant for the 
photo- and electroproduction of charmonium mesons are depicted in Fig. 4. There, we 



employed nonrelativistic charmonium wave functions calculated from a power-law p9[ and 
a logarithmic potential [0, which both describe Fj/^^e+e- reasonably well, and we set 

-, with /i = 



1 + ^ ) . This yields 



(*)- 



r oc 



ml + for the wave function 



of the cc component of the (virtual) photon. Note, however, that at sufficiently large 
- which is the only domain at which the formulae for any vector meson production can 
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be consistently proven in QCD - / (Pktipj/^{kt, z) oc z{l — z) due to the Q"^ evolution 
5I|,[TP|, which is qualitatively different from the behavior expected for the nonrelativistic 



charmonium wave functions. The question which arises is thus whether hard physics is 
applicable at smaller Q^, and so we are effectively analyzing physics hidden under the 
assumption that the nonrelativistic charmonium wave functions approximate the light- 
cone wave function. 
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FIG. 4. Average transverse sizes of the cc components effective in ^^(*)Ar^viAr for J/^* and 
photo- and electroproduction. 



We observe that for J/^ photoproduction the corresponding relevant distances, h^, 
are similar to the average interquark distances in the J/\& mesons (~ 0.4 fm), while for 
the photoproduction of the excited state, the \E'', the distances that are most relevant 
for the photoproduction process are significantly smaller than the respective average in- 
terquark separations (~ 0.7 fm). Within the charmonium model approximation for the 
light-cone wave function of the \t'' meson, this is due to the nodal structure of the wave 
function of the \E'', as was also observed by Benhar et al. ||l6l in the framework of the 



two-gluon-exchange model |18|. However, for charmonium electroproduction at larger Q 
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the distances relevant in the amphtude are much smaller than the average interquark 
separations in the charmonium mesons. 

We can now also estimate the contribution of the hard blobe to the t dependence of 
the vector meson electroproduction cross section. In detail, we fit 

(fb dz i)^, (z, b) a{b^) ^viz, b) e-^"P*^ oc e^^*/' _ ^ (20) 



by an exponential at the vicinity of t = 0. The slope parameter that we find, By ^0.7 
GeV~^, is much smaller than, e.g., the experimentally observed slope. Bp = 4.6±0.8 GeV~^ 



3^ , of the p-production cross section off the nucleon. Hence, for all current practical 



considerations, the hard blobe does not contribute significantly to the t dependence of 
the cross section, which hence originates almost entirely from the target's two-gluon form 
factor, G2g{t), and which is thus expected to be universal for all hard processes driven by 
the gluon density. 

In the case of a transversely polarized 7*, the contribution of large b, where nonper- 
turbative QCD effects are dominant, is not suppressed. Those contributions arise from 
very asymmetric qq pairs with z ~ or 1 p3|,p]JT^, and, for transverse polarizations, in 
the integrand of the amplitude the regions z or (l — z) ^ 1/Qb thus give important contri- 
butions, which are usually suppressed at large by a Sudakov type form factor. In fact, 
it is easy to check that, for the case of transversely polarized photons, the contribution 
of this regions to the asymptotic expressions contains a logarithmic divergency. Hence, in 
this case, the integral over transverse momenta does not decouple. 

On the other hand, the increase of xGn{x, b) at small x works towards enhancing the 
perturbative contribution of small b in this case as well. This is because the contribution of 
quark configurations with small b increases faster with decreasing x than the contribution 
from soft QCD. At what x and the pQCD regime will dominate for reactions initiated 
by a 7^ is an open question since the contributions from soft QCD are not under control in 
this case. Currently, we are nethertheless estimating the hard, perturbative contribution 
to a^^N-^vN 131. 

Furthermore, since the formulae in this paper were deduced in the leading as In -r^ — 
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approximation only, there is an uncertainty in the theoretical formulae whether the gluon 
distributions should be evaluated at or say at (5^/2. A similar problem exists for 
calculations of the leading twist contribution to the total DIS cross section within the 
leading In -r^ — approximation. Note, also, that the substantial difference between 
ba^iQ"^) and b^^^viQ^) which we found in this subsection indicates that significant next- 
to-leading order corrections should be present in diffractive vector meson photo- and 
electroproduction. We will give phenomenological estimates of these effects in Sects. IV. C 
and V. 

B. The nondiagonal parton distributions 

In this subsection, we show that the nondiagonal parton distributions of the target can 
be expressed through the diagonal - conventional - parton densities. The first observation 
is that the momentum transfered to the target, tmin = —m'j^x'^{l + is negligible. 

A more important issue is that the difference between the masses of the virtual photon 
and the vector meson leads to a difference between the light-cone fractions of the target 
momenta in the initial (xj) and final (x/) states. To elucidate this point, let us consider 
diffractive vector meson electroproduction in the center of mass of the 7* -|- T system. 
Suppose the initial target has momentum P. Kinematical considerations show that the 
final target has momentum P(l — {Q^ + My)/s), and thus the difference between Xi and 
Xf is X p. 

For the hard part of the diffractive amplitude, the dependence on fractions of the target 
momentum carried by the interacting partons in the initial and final states is calculable 
in pQCD. This simplification arises since in the kinematics of diffractive processes, which 
dominate in the considered diagrams, it is legitimate to neglect in the propagators of 
the exchanged particles the contribution of the longitudinal part of their momenta as 
compared to the contribution of the transverse components. This can be easily proven 
within the agln-r^ — and/or In -r^ — In - approximation by direct evaluation of the 

^^QCD ^^QCD ^ 

leading Feynman diagrams. In reality, this argument is valid beyond the pQCD regime 

19 



since it explores specifics of multiplieriplieral and multiregge kinematics. 

Tlius, we conclude that the initial and final parton wave functions of the target practi- 
cally coincide, and it is a reasonable numerical approximation to use the diagonal parton 



densities. This question will be considered in detail in Ref. ||25|. Here, we have explained 
that the tiny momentum transfered to the target does not influence the hard blobe, which 
is specific for small x physics. 



IV. ESTIMATES OF PREASYMPTOTIC EFFECTS 

A. Vector mesons built of light quarks 

The calculation discussed above is applicable at suffiently large Q^, i.e., only when the 
transverse momenta in the qq pair are neglected, can the cross section for electroproduction 
of longitudinally polarized vector mesons be written in the form of Eq. (^). To analyze 
the onset of the asymptotic regime, we consider here the effects of the Fermi motion of 
the quarks within the diffractively produced vector meson. As originally suggested in Ref. 
p5| , we do this by keeping the quarks' momenta in the propagator of the virtual photon, 
as the numerical coefficient which accompanies the leading term is large (~ -^)- We 
will show that account of Fermi motion of the quarks within the vector meson leads to a 
significant suppression of the cross section at moderately large Q^. Thus, the formulae for 
diffractive electroproduction of vector mesons accounting for this Fermi motion contain 
an additional factor, T^Q"^), which can be straightforwardly obtained from Eq. (2.15) of 
Ref. [ffl, and which, within the approximations discussed above, has the form: 



V 



(21) 



Here, At is a two dimensional Laplace operator which acts on the wave function of the 
photon. We neglect here the current quark mass of the light quarks, which is small on the 
scale of the discussed phenomena. 
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To estimate the role of Fermi motion effects, we made estimates of the dependence 
of T for three types of vector meson wave functions, ipviz, h) ~ z{l — z)(f){kt). For the 
nonperturbative piece, (f){kt), we choose either (l){kt) oc jp^jjTj2, as in the last section (la- 
beled "dip."), or0(fc t) (X exp{—kt/ fi), as motivated by the observed transverse momentum 
distribution of direct pions in hadroproduction (labeled "exp."). We consider also a model 
where (j){kt > IGeV/c) oc (labeled "reg.") in order to check the sensitivity of our re- 
sults on the large transverse momentum tail of the nonperturbative wave function. Such 



a behavior (without the naively expected as{k^)) follows directly from QCD |36|. At the 
same time, the region of applicability of the pQCD prediction for the high momentum tail, 
which we started at 1 GeV/c, is still subject of discussions. We then vary fi in the range 



corresponding to average transverse momenta y ~ 300 -i- 600 MeV/c. The results of 
these calculations are presented in Fig. 5. 

1.0 
0.8 
0.6 



300 MeV (exp.) 
300 MeV (dip.) 
300 MeV (reg.) 
600 MeV (exp.) 
600 MeV (dip.) 
600 MeV (reg.) 




0.4 
0.2 h 
0.0 

1 2 10 20 100 

[GeV^] 

FIG. 5. Suppression factor, T(Q^) of Eq. ([2l]), in daj*^j\[-,vN /dt\t=o due to transverse quark 
Fermi motion for diffractively produced vector mesons built of light quarks. 



One can see that a very signicant suppression of the vector meson yield may be present 
in the range covered so far by the HERA experiments, and that the high momentum 
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tail - due to hard interactions between the constituents of the vector meson - seems to be 
an important correction in the relevant range of ^ 10 GeV^. As can be seen from Fig. 
5, at large Q^, T is determined mostly by the high momentum tail of the kf distribution 
of the bare quarks within the diffractively produced vector meson. In mean field models 
of hadrons - bag models, constituent quark models with an oscillator potential or models 
based on the method of dispersion sum rules - this tail is included into the effective 
parameters, and thus these models give no clues for the value of this component, or its kt 
and z dependence. On the other hand, as suggested by Fig. 5, a study of vector meson 
production over a wide range of would allow to obtain unique information on the kt 
distribution in hadrons. Note, however, that to achieve a nonambigous interpretation of 
quark Fermi motion effects, it is necessary to calculate also the contribution of the qqg 
component in the light-cone wave functions of the 7^ and the vector meson to the cross 
section of hard diffractive processes. This has not been done so far. 



In the case of the photo- and electroproduction of heavy vector mesons, the overlap 
integral is suppressed both due to the presence of a term proportional to ml and due to 
quark Fermi motion. 

To investigate the role of Fermi motion for a realistic case, we consider the non- 
relativistic model of charmonium. We keep terms proportional to ^ in the wave 
function of the photon, as they enter with a large numerical coefficient and because 
the integral over the charmonium wave function only slowly converges at large mo- 
menta. At the same time, we neglect relativistic effects in the charmonium wave func- 



^This formula follows from the necessity to transform light-cone perturbation theory diagrams 
into conventional, nonrelativistic diagrams describing the charmonium bound state. 



B. Charmonium production 




In this approximation, the phase volume in, for 
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instance, Eq. 



jT|) becomes 



into 



=d^k, the energy denominator transforms 



^(1-2) 



4(m^ + k^), and the charmonium's hght-cone wave function %l)y{z,kt] 



can be obtained from its non-relativistic rest frame wave function ipvik) based on the 
requirement that the J/\l/ — e^e~ decay width is reproduced correctly, which yields 
/ dz d'^kfipviz, kf) = J d^kipv{k)- In this limit, rjy ^ 2.44, and the corresponding factor, 
Fy, is significantly smaller than unity: 



z(l-z) 



V 



(22) 



Here, At is the two dimensional Laplace operator which acts on the wave function of 
the photon. For numerical estimates, we use two different realistic charmonium wave 
functions calculated from a QCD-motivated ||3^ and a logarithmic potential 0], which 
both describe Tj/qi^e+e- reasonably well. We hereby restrict our consideration to potential 
models for which the mass of the constituent c quark is close to the mass of the bare current 
c quark, i.e. rric ^1.5 GeV {rric = 1.48 GeV for the QCD-motivated potential of Ref. |P7|| , 



and rric = 1.5 GeV for the logarithmic potential of Ref. 0]). This is necessary to keep a 
minimal correspondence with the QCD formulae for hard diffractive processes which are 
expressed through the distribution of bare quarks [Q . The results of these calculations are 
shown in Fig. 6 for diffractive electroproduction of J/\l/ and mesons, respectively. Fig. 
6 shows that, in this model, the applicability of the asymptotic formulae requires rather 
large Q'^. 

Note that, in the evaluation of rj/,j,^e+e- and also in the process considered here, 
large radiative corrections, 1 — ^^"3^"^"'' , accompany the charmonium's non-relativistic 
wave function but not the qq component of the light-cone wave function, which underlines 
the limits of applicability of non-relativistic charmonium potential models in that context. 
Inclusion of these radiative corrections effectively corresponds to a renormalization of the 
nonrelativistic wave function by the factor 1— ^""^^"^ Since we only employ the light-cone 
normalization of the wave function, this factor is implicitly included in all our calculations. 
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FIG. 6. Suppression factor, FyiQ"^) of Eq. (|2^), in da^^N-^vN /dt\t=o due to the finite quark 
mass and quark Fermi motion for diffractively produced J/^ and ^I*' mesons. 

Furthermore, the nonrelativistic charmonium model does not include gluon emission at 
higher resolution, and it therefore does not display the expected asymptotic behavior: 
/ (Pktipviz-, kt) oc z{l — z). In principle, this could be accounted for by means of 6-space 
evolution of the leading twist wave function as suggested by Brodsky and Lepage 

In Ref. pQCD was applied to the diffractive production of J/'^ mesons starting 
from photoproduction. The observed increase of the J/\E' photoproduction cross section 
with initial energy is in line with pQCD expectations . Moreover, the theoretical analysis 
of the transverse distances between the c quarks within the J/'^ meson shows that they 
are close to that dominant in the wave function of the 7* in the deep inelastic structure 
functions. Thus, the necessary condition for the applicability of pQCD is the same in all 
cases. Specific to charmonium, however, is that these distances are close to the average 
distances between the quarks within the J/'^ meson, as was discussed in the previous 
section. So, the process under consideration should be very sensitive to the exact form of 
the light-cone wave function of the J /'^ meson and to its color field. 



J/V- [QCD-motiv.] 
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It was assumed in Ref. that the cross section of diffractive J/\I/ electroproduction 
could be calculated employing the nonrelativistic constituent quark model for the J/\E' 
meson's wave function, both for transverse and longitudinal photon polarizations, and 
quark Fermi motion was neglected. This corresponds to assuming that R = rij/ii,Tj/ifi = 2. 
In the limit where we can justify the application of pQCD (Mj^^ ^ Q^), our result 
coincides with that of Ref. if we neglect Fermi motion and 7^ Mj/,j,/2 effects, i.e., 
assume that R = 2, and if we discard the qualitative difference between the minimal 
cc Fock component in the J/\E''s wave function and its wave function in nonrelativistic 
potential models, which is actually not legitimate as was discussed in Sect. III. A. Thus, 
in Ref. 0, the predicted cross section is overestimated within the nonrelativistic quark 
model by a factor Ty^{Q^), where 

2 



TviQ' 



(Q2+4mg)2 



V 



(23) 



which is presented in Fig. 7 for the charmonium wave functions employed previously. 

One can see from Fig. 7, that the effects of Fermi motion strongly suppress the 
contribution of the cross section, especially for ~ where most of the experimental data 
were obtained. If we evaluate the suppression factor Ty of Eq. ( p3D at the photoproduction 
point with charmonium wave functions calculated from the two nonrelativistic potentials 
of Refs. 0,^], we find a value of Ty(Q^ = 0) ~ 1/9. This is significantly smaller than 
the corresponding value (~ 0.5) given recently by Ryskin et al. 



In our opinion, the reason for this discrepancy is threefold: Firstly, the authors of Ref. |38| 



use the approximate Fermi motion correction which we gave in Ref. |35|, secondly they neglect 
the longitudinal relative motion of the quarks in the J/^' mesons - by setting z = ^ - which, 
in turn, transforms the three dimensional integral over d^k in Eq. (^) into a two dimensional 
integral over d'^kt - see Eq. (23) in Ref. [p8| - and, thirdly, a Gaussian form for the wave 
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FIG. 7. Correction factor, Tv{Q^) of Eq. (p3|), in da^^N^VN /dt\t=Q for diffractive charmo- 
nium production due to quark Fermi motion. 



Tv{0) has been evaluated also by decomposition over powers of Ty(0) = 1 



20 / \ 
3 \m2/- 

is due to 



The difference between the actual number of 20/3 and the 4 obtained in Ref. 
account of all terms arising from the application of the operator. By definition, {k"^) 



Ji>vik)d3k ■ 



In Ref. Gaussian type wave functions for the J/ip were used to estimate 



that (k'^/ml) ^1/8 (Note that for a pure Coulomb potential, which, when masked with 
a running coupling constant, corresponds to the limit rric — ^ cxo in QCD, (fc^) = oo due to 



functions was assumed in Ref. [^8[, whereas the realistic potentials of Refs. [^,3C,37] all yield 
wave functions that fall off significantly more slowly. All three approximations diminish the 
role of large quark momenta and result in an essential reduction of the suppression factor. For 
instance, our TviQ"^ = 0) would rise from 0.11 to 0.4 if we were to use our approximate formula 



of Ref. [35 1 and change the integration over d^k in Eq. ( P3|) to an integration over d'^kt, and 
it would rise further to 0.5 if we would replace our realistic charmonium wave functions with 



Gaussians. Note that the latter value coincides with that given in Ref. [38|. 
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a linear divergency at large k.). So, within this approximation, T'y(O) ~ 1/(1 + a^)^- 

The suppression factor TylQ"^) has contributions from leading and non- leading twist. 
The respective leading twist expression follows from Eq. (pTf ) by setting, in the wave 
function of the photon, kt to zero after differentiation: 



^ ^ z{l-z) . 



..2 



V 



(24) 



Additional (leading twist) terms containing Inm^ will be accounted for in a separate 
publication. Our analysis shows that, in QCD, the higher twist corrections are smaller, 
which, in turn, justifies our discussion of those Fermi motion effects in what is otherwise 
strictly a leading twist analysis. 

Furthermore, the results of calculations of J/"^ and photoproduction cross sections 
are rather sensitive to the value used for the mass of the c quark, and, for instance, in 



Ref. 1^ the assumption was made that 



rrir. 



M J 1^/2. If we analyze the dependence of 



the J/\l/ photoproduction cross section on the difference between Mj/,j, and 2mc, we find 
a correction factor, 

2 

(25) 



R 



M 



( 



2mr 



Jd^k 



Jd^k- 



k'^+rrA 



(l+4A:2/Mj/^)2 k'i+ml 

which, using the values for iric and the respective charmonium wave functions from Refs. 
5I1| , |57[ ], yields RM,/^j^2mc ~ 1-2- This, together with our Tv{Q^ = 0), gives a final predic- 



tion for the suppression of J/\l/ photoproduction olQ ~ 1/8. 



^In Ref. 1^^, the correction factor of Eq. (^) was approximated by {Mj^q,/2mc)^ , and the 
latter was evaluated with a small rric = 1-43 GeV, as estimated by Jung et al. ||3^, which, in 
turn, yields a very strong enhancement of RMj^^^2mc ~ 2. This estimate is at variance with 
the direct calculation using Eq. (pst). This, together with the underestimated Fermi motion 
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C. rescaling 



In Sect. III. A, we determined the relevant transverse distances between the quark and 
the antiquark in the hard blobe of Fig. 1, and we found that they depend on the process 
under consideration, i.e. there is a substantial difference between b^^{Q^) and b^^^viQ"^)- 

In detail, we used the impact parameter representation of the longitudinal structure 
function, aii^x, Q^) of Eq. (|18D, to define a relation between the transverse size of the qq 
pair, b, and the virtuality of the process, Q^. This relation. 



b = b.AQ')^^^^ (26) 



at a; = 10^'^, yields the scale effective in the gluon densities which enter, for instance, 
in the elementary color-dipole cross section, of Eq. (|^. It is designed in such a 

way that, if we were to evaluate the integrand of Eq. (^) at this average b = 60-^ (Q^) 
instead of carrying out the integral, we would just recover the normalization condition: 
<7l{x,Q'^) oc as{Q'^) xGn{x,Q'^). Hence, the virtuality that corresponds to the dominant 
transverse distances and the virtuality of the process are identical, at least for the case of 
the longitudinal photoabsorption cross section. 

However, this is not the case for the electroproduction amphtude, A^^n^vn of Eq. 
(P!B|), due to the significant difference between the transverse distances it is dominated by 
and bcr^iQ"^)- This indicates that substantial next-to-leading order corrections should be 
present in diffractive vector meson photo- and electroproduction which, in turn, should 
lead to a change of the scale effective in the gluon densities, i.e., the virtualities that enter 
in the argument of as^Q"^) xGn{x, Q"^) in our basic equation (|T^) . 

Following our logic of approximating the integrals over the transverse sizes of the qq 
pair effective in the hard blobe of Fig. 1 by the contribution from an appropriately defined 



correction, leads the authors of Ref. |3Sl to the conclusion that there is no net suppression for 
J/^* photoproduction. 
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average transverse distance, a natural guess is to rescale the virtuality in proportion 
of these average transverse sizes, 

^'^^ ^^ZW) ' ^^^^ 

which, for p production, yields significantly smaller virtualities, i.e., Qe// < Q^- This 
results in a strong decrease of the quantity [a;s(Q^) xG'7v(a;, Q^)]^, and therefore also of 
the p-meson electroproduction cross section, as can be seen in Fig. 8 where we show 
\as{Q'^) xG]^{x,Q'^y^ with and without the NLO rescaling.[] 

In Fig. 8, we also plot the exponent of the energy dependence of [as(Q^) xGj^{x, Q^)]^ 
(and hence of the predicted cross sections), ma;(Q^), which is defined by means of 
[as(Q^) xG'Ar(x, Q^)]^ oc x~™"''^'^'^\ We observe that the rescaling not only significantly 
reduces the predicted cross sections, but it also results in a slight reduction of its steep 
energy dependence. 

On the other hand, the relevant transverse distances for J/\E' photoproduction, 
= 0), are relatively small, = 0) ~ 0.25 fm, which, using the logic outlined 

in the above, results in substantially higher effective virtualities, cf' ^ ~ 5.1 GeV^, 
for J/\l' photoproduction than those suggested, for instance, in Ref. - Q =2.4 GeV^. 
This will be discussed in more detail in Sect. V.C 

In Fig. 9, we finally show the Qe// obtained in this manner for the cases of p-meson 
electroproduction and J/\l/ photo- and electroproduction. Note that the relationship be- 
tween and Qe// is almost independent of x, and that the Qe// production 

differs significantly from the estimate of Ref. Q = which is shown through 

the dotted line in Fig. 9b. 



Eq. (27) assumes that Gn{x^Q ) is defined such that Eq. ^ for is valid in the NLO 
approximation. This question requires further investigation. However, the corresponding un- 
certainty is numerically unimportant for the comparison of production rates of different vector 
mesons. 
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FIG. 8. The square of the gluon density, [asiQ"^) xGn{x,Q'^)\^ , and its energy dependence, 
m^{Q^), defined by [asiQ"^) xGN{x,Q'^)f oc calculated for the CTEQ3L parton distri- 

bution function with and without the NLO (/O-meson electroproduction) rescaling introduced in 
Eq. d^). 
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FIG. 9. Qlff obtained via the NLO rescaling introduced in Eq. ( p7|) for p-meson electro- 
production and J/^* photo- and electroproduction. For the latter, we also show the estimate of 
Ref. pi (the dotted line). 
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V. COMPARISON WITH DATA 



A. Main features of the data 



The HERA data 0, when combined with the NMC data on diffractive electro- 
production of p-mesons, are consistent with several key predictions of Eq. (p!^): 

(i) a fast increase with energy of the cross section for electroproduction of vector 
mesons, which is approximately proportional to x~^'^ for = 10 GeV^. Note that this 
fast increase with decreasing x is not expected within the non-perturbative two-gluon- 
exchange model of Donnachie and Landshoff |^ , which predicts a cross section rising as 
~ 3;-o.i6 ^ _ even weaker increase of the cross section when integrated over 
t due to a substantial increase of the slope of the elastic cross section expected in soft 
Pomeron models: aei Qc ;;;cT6 gp^2a°iii^ " Here, Bq is the slope of the elastic cross section 
for X = Xq. 

(ii) the dominance of the longitudinal polarization, oc P,^. Data do confirm 
the dominance of ai, though they are not accurate enough to study details of the 
dependence. Predictions for the x-dependence of this ratio sensitively depend on the in- 
terplay of perturbative and nonperturbative contributions to ax- In the leading as In -rS — 

^QCD 

approximation, it is easy to generalize the QCD evolution equation to transversely po- 
larized photons as a projectile As in the case of the regular structure functions, the 
appropriate question will be on the region of integration over fractions of momenta appro- 
priate to the hard blob. Since the hard kernel in the evolution equation increases faster 
at small x than the soft contribution, one may expect an onset of hard physics for trans- 
versely polarized photon projectiles as well, though at smaller x and at higher Q^. Note 
that at sufficiently large Q^, the contribution of nonperturbative QCD to diffractive vector 
meson production is suppressed in addition by a Sudakov type form factor. At sufficiently 
large Q^, the process should be discussed in terms of the Bethe-Salpether wave function 
for the transition of a p-meson into two highly virtual quarks. Such a mathematical object 
exists in pQCD only. Thus, there exists a possibility to check quark counting rules in their 
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original form. This may manifest itself in the x-dependence of the ratio ^ for fixed Q^. 
At intermediate ~ 10 GeV^, where hard physics already dominates in ctl, gt may 
still be determined by soft, nonperturbative contributions. For these Q^, the ratio should 
increase with decreasing x as x^G\[x^ Q^). At high Q^, where hard physics dominates for 
both Gi^ and cr^, the ratio would not depend on x. We are currently investigating these 
issues in detail [Q. The x dependence of ^ was also discussed in Ref. The authors 
employed nonrelativistic oscillator wave functions for the diffractively produced vector 
mesons, and the pQCD color-dipole picture was applied for both polarizations and also 
for large transverse distances. For both cases, the formula for the elementary color-dipole 
cross section, cr(&^) of Eq. (|^) which was deduced in |21|J2^, is clearly not applicable, as 
was discussed in detail in Sect. III. 

(iii) the absolute magnitude of the cross section is reasonably well explained, though 
uncertainites in the knowledge of xGAr(x, Q^) are rather large to make a comparison 
quantitative. This will be discussed in detail in the next subsection. 

(iv) the dependence observed by the ZEUS collaboration: ctl oc with n = 
4.2 ± 0.8 0]. This seems to reflect the evolution of the parton distributions as well as 
effects of quark Fermi motion within the vector mesons, as will be outlined in subsection 
V.D. 



B. Absolute magnitude of the p-electroproduction cross section 

Currently, absolute cross sections for exclusive p-meson production are available from 
NMC m and from ZEUS [§. 

We restrict our analysis here to > 6 GeV^. In the case of NMC, the most relevant 
data are those obtained using a deuteron target as no data were taken with hydrogen. 
The measurements do not separate the elastic deuteron final state and deuteron break up. 
So, to extract the cross section for p-production off the nucleon from the deuteron data, 
we use the closure approximation to sum over the final states of the two nucleons, which 
yields [P 
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= 2(1 + F,,(4t))^^^^ , (28) 

where F£,{t) is the deuteron's electromagnetic form factor. In the above, we neglected 
shadowing, Glauber rescattering effects (which are suppressed due to color transparency 
expected in this kinematics) and the real part of the elementary amplitude. Fu(t) can be 
roughly parameterized as Fo{t) ~ c^-^* with Bo ^ 14 GeV~^ |^T|. This yields 



5D 

where B = 4.6 ± 0.8 GeV^^ [^] is the slope of the p-production cross section off the 
nucleon. With this, the total cross section of p-production on the nucleon, a^*N^pN, can 
be calculated from the NMC's deuteron data. Using their determination of the virtual 
photon's polarization, (e), and the measurement of i? = 0"^.7v^pAf/'^7*7V^pAf = 2.0 ± 0.3 at 
(Q^) = 6 GeV^, we can furthermore separately extract the elementary longitudinal cross 
section, ai:j*^_^p^. In the same manner, a^^^^p^ can be obtained from the ZEUS data 
where i? = 1.5 ± 0.6 (statistical errors only) at {Q"^) = 11 GeV^. 

Results of such calculations are presented in Fig. 10. There, we also show theoretical 
predictions of Eq. ([T^ employing several recent parametrizations of the gluon density]^ 
^,43,^ and typical parameters for the p- meson wave functions as considered in section 



IV. We set r]y = ?) and parametrize the dependence of the differential cross section on 
the momentum transfer in exponential form withj^ B A GeV~^. Note that a change of 
T{Q'^) in the range depicted in Fig. 5 introduces an extra scale uncertainty of 0.5 ^ 2. 



^Note that, within the leading In — approximation to QCD, the use of NLO fits to the 
nucleon's gluon density is not fully consistent. 



natural estimate for B is B = -^^^ + B^*^p{Q'^) where \J {^1) is the mean radius of the gluon 
distribution in a nucleon which is usually taken in realistic models of the nucleon as 0.6 fm. The 
calculation of Sect. III.A gives B^^p^Q"^) ~ 0.7 GeV~^. Putting all numbers together we obatin 
B w 3.8 GeV-2. 
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FIG. 10. The total longitudinal cross section, ai^tj^^pj^, calculated from Eq. (|T2|) for several 
recent parametrizations of the gluon density in comparison with experimental data from ZEUS 
|] (full circles) and NMC |3| (squares). 
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We discussed in Sect. IV. C that next-to- leading order corrections should lead to a 
decrease of the scale effective in the gluon densities, i.e., the virtualities that enter in the 
argument of as{Q'^)GN{x,Q'^), due to the fact that 6^£^p((5^) is substantially larger than 
^cri(Q^)- This results in a decrease of the p production cross section as illustrated in Fig. 
10 for the CTEQ3L parametrization through the dotted lines. 

One can see from the figure that we substantially underestimate the cross section at 
lower energies which, in turn, correspond to rather large x: for example (x) = 0.103 for 
the NMC point at = 16.9 GeV^. This kinematics is well beyond the x-range specified 
in Eq. (H), and the frozen approximation, which we use, is inapplicable. Qualitatively, 
one can expect that, in this situation, quantum diffusion of the qq pair would lead to an 
increase of its interaction cross section, and hence to agreement with the observed cross 
section. 



C. Absolute magnitude of the J/^' photoproduction cross section 



In a recent work by Ryskin et al. |^8[, J/^ photoproduction data were compared 
with the corresponding pQCD predictions. For the latter, the nonrelativistic constituent 
quark model for the J/\E''s wave function was employed, and Fermi motion and Mj/^ ^ 
2mc corrections were estimated to leave the cross section unaltered. On the other hand, 
we argued in Sect. IV. B that there is a strong reduction (~ 0.13) of the naive pQCD 
prediction for J/\E' photoproduction due to those presymptotic effects. Here, we thus 
compare our predictions with the data, while either accounting for the preasymptotic 
corrections, i.e. the factor Ty(Q^) of Eq. (^) and the Mj/^ ^ Irric correction of Eq. 
or while leaving them out. 



Setting, r]y = 2 and 



^2 



Q 



as suggested in Ref. our Eq. (|T2D yields 



dt 



-Mj,Ws{Q)Tv{Q^ 



t=o 



d 



2 dlnx 



(30) 



12aEMQ 

which agrees with the result given in PJ3^ if we neglect the Fermi motion and Mj/^ ^ 2mc 
corrections, i.e., if we set TyijQi'^) of Eq. (|23| ) to 1 and replace rric by Mj/^/2. On the 
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other hand, we can also use the results of our analysis on the relevant transverse distances 
of Sect. III. A (see Fig. 4) to estimate the scale, Q, effective in as{Q)xGN{x,Q) in Eq. 
(^). As outlined in Sect. IV. C, we find for J/\E' photoproduction Q ^ 5.1 GeV^ which 
is substantially higher than the estimate of Ref. ||38| - Q = 2.4 GeV^ - and which leads 
to an increase of the respective cross sections by a factor of about 1.5 for the logarithmic 
potential of Ref. [Q. Taking the experimental value for the slope parameter of -Bj/* ~ 4.5 



GeV~^, the cross section for diffractive J/\l/ photoproduction is calculated from Eq. (|30|), 
and corresponding results are shown in Fig. 11 for various recent parametrizations of the 
gluon density. 

Thus, if Fermi motion, the 7^ correction and rescaling are accounted 

for, the predictions of the pQCD model and the photoproduction data agree, at large W, 
within the uncertainties in the existing gluon distributions for the logarithmic potential 
model pO|, which gives a smaller suppression as it has a rric close to the current quark 
mass (cf. discussion in Sect. IV. B). 

It is instructive, also, to consider the ratio of the J/\& photoproduction and the longitu- 
dinal p-meson electroproduction cross sections, cr(7j9— * J/\l/p)/cr(7£p— >pp), for = 8.8 
GeV^ and at the same x. We find that the various gluon distributions lead to very sim- 
ilar ratios since the gluon densities enter at about the same virtualities. In detail, we 
find cr(7p^ J/\l/j9)/cr(7£p— s>pj9)|^;jgQ ~ 0.9 with the Q rescaling and 0.6 without both at 
X = 2.1 ■ 10^^ and x = 8.5 ■ 10^^ with a theoretical error ^ ±0.4 due to uncertainties of 
the estimate of the suppression factor Ty(Q^). The experimental data yield 1.2 ± 0.6 at 
X = 2.1 ■ 10"'^ and 1.6 ± 0.7 at x = 8.5 ■ 10~^, and thus the theoretical prediction for this 
ratio is no more than la below the value we extract from the data. 

Note, also, that, in a truely relativistic description, for photoproduction - where only 
transversely polarized photons contribute - the suppression factor of Eq. (p3D would 
contain an extra factor oc 4^(Y^' "W'hich, in turn, would reduce the suppression by about 
20%. At the same time, the analysis of the essential interquark distances (in Sect. III. A) 
seems to indicate that the coupling of the two gluons to the charmonium meson is more 
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FIG. 11. The diffractive J j'^ photoproduction cross section compared with experimental 
data from E401 [|l^ (full circles), E516 |ll| (diamonds), ZEUS '93 (squares), and preliminary 
ZEUS '94 (circles) and HI (stars) data |l|l. (a): no corrections (Q = Mj/^/2), (b): with the 
corrections of Sect. IV. B (Q = w-c) (c): ditto (but with rescaling: OsxGn evaluated at Q 
based on bj^^). gg 



complicated than assumed in the nonrelativistic description. In particular, the interaction 
of the exchanged gluons with the interquark potential is not suppressed. Only in the 
approximation where the mass of c quark tends to infinity, this effect would be small 
(~ where It is the transverse momentum of the exchanged gluon. 



One should also notice that in the fixed target experiments 0,|TT|], the recoil proton 
was not detected and the corresponding cross sections could thus have been overestimated. 
Also, it was demonstrated by Jung et al. ||39| that, in the kinematics of these experiments, 
the photon-gluon fusion mechanism could yield up to 50% of the cross section. Besides, 
the coherence condition of Eq. is violated for J/^ photoproduction for W < 12 GeV. 

Also, the qualitative agreement of our predictions for the J/\l' photoproduction cross 
section with the data hints that the same approach could be used to calculate the total 
J/\E'-N interaction cross section. The inspection of Fig. 3 - at 6 ~ 0.40 fm - shows that 
such a model approach predicts a rather fast increase of this cross section with energy (in 
this case x = Mj^^^f s) which may appear important for the interpretation of AA collisions. 
Note that the transverse distances relevant for J/\E' scattering are larger than those for 
J/^ photoproduction where small 6's are enhanced due to the presence of the photon's 
wave function in the integral. 



D. dependence of p-meson electroproduction 

Let us now analyze the predictions of Eq. (|T^ for the dependence of the vector 
meson yield. One may conclude from inspection of this equation that it leads to a 
dependence of the cross section. There are, however, two QCD effects which substan- 
tially modify this naive expectation. At large Q^, where higher order corrections can be 
neglected, Eq. (|12]) predicts a dependence of the cross section which is substantially 
slower than because the gluon densities fastly increase with at small x. Numer- 
ically, the factor {as{Q) xGn{x,Q))'^ in Eq. (|T2|) is oc with ~ 1. To illustrate this 
point, we define the effective power of the dependence of [as{Q) xG n{x, Q)]"^ as 
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meff{x,Q'^) = 2 



d\n{as{Q)xGN{x,Q)) 
d\nQ 



(31) 



Results for meff{x,Q'^) using several current parametrizations of the gluon distribu- 
tions in the nucleon |]42| , ^ , |27[| are presented in Fig. 12 for (5^=10 and 40 GeV^ as functions 
of a; (dotted lines) together with the experimental ZEUS data point 0. One can see that, 
typically, they reduce the dependence by one power of Q for Q"^ ~ 10 GeV^. 

A preasymptotic effect may come from higher order corrections to the factor, see, 
for instance, the discussion in Sect. IV. The corresponding term, T{Q'^) of Eq. ( pTD and 
depicted in Fig. 5, yields an effective power of peffiQ"^) = ^^gfj^q ^ to the dependence 
of the cross section, where pe//(10GeV^) ~ 1.5 and pe//(40GeV^) ~ 0.8. Thus, combined, 
these two effects result in 



as shown also in Fig. 12. One can see from the figure that, in the ZEUS range, the rie// 
we find is consistent with the experimental number of ZEUS. Note that if we include 
rescaling, which was discussed in Sect. IV. C, the ^e//'s we find increase by about 0.3 for 
the kinematics of the HERA experiment. 



Longitudinal vector meson production is dominated by small interquark distances in 
the vector meson wave function. Therefore the factorization/decoupling theorem can be 
used to calculate the cross section for hard diffractive processes in QCD without any 
model assumptions. For My <^ Q^, all dependence on the quark masses and thus on 
flavor is contained in the light-cone wave function of the vector meson and not in the 
scattering amplitude. This prediction is non-trivial, since experimentally the coherent 
photopro duct ion of mesons containing strange or charm quarks is strongly suppressed as 
compared to the SU(4) prediction for the ratio of the production cross sections for various 
vector mesons, which is 




(32) 



VI. RESTORATION OF FLAVOR SYMMETRY 
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p° : uj : 4) : J/^ = 9:1:2:8 



(33) 



Experimentally, this suppression factor is ~ 4 for 0-mesons and ~ 25 for the J/\&. Thus, 
QCD predicts a dramatic increase of the : p° and J/\l/ : p° ratios at large 0. We 
expect that for the (p : p-ratio, the pQCD limit would be reached rather early, while the 
restoration of SU{A) in J/\E'-meson production would require extremly high Q^. This 
will be clarified further in the next section. Besides, QCD predicts a slow increase of 
the relative yield of heavy flavor production with decreasing x since the virtuality which 
enters in the gluon distribution is larger for heavy flavor production. 

Moreover, an additional enhancement of heavy flavor production is expected since, for 
heavy quarkonium states, the probability for a q and q to be close together is larger. In 
fact, Eq. (|T^), which was derived from QCD, predicts for the production ratio of mesons 
Vi and V2, at large Q^, that 

a(72 + r-. Vi+T) 



aiYL + T^V^ + T) 



My, Tv.^e+e- VW) (34) 



Values of rjv for mesons built of light quarks seem to be close to the asymptotic value of 
rjy = ?) already for moderate ~ 10 GeV^. This reflects the observation of a QCD sum 
rule analysis as well as lattice QCD [^] and the instanton liquid model that the 
nonperturbative interaction in the vector channel is weak. 

Based on the measured values of Vy^e+e- and estimates of rjy for p and and the 
charmonium model for J/\I^, we observe that Eq. (|l^) predicts a significant enhancement 
of and J/\l/ production as compared to the SU(4) prediction of Eq. (|33|): 

p° : : : J/^ = 9 : (1 * 0.8) : (2 * 1.0) : (8 * 1.9) , (35) 

At very large 3> My, the qq wave functions of all mesons should converge to a universal 
asymptotic wave function with rjy = "i [^]. In this limit, further enhancement of heavy 



resonance production is expected, 

p° : : : J/^ = 9 : (1 * 0.8) : (2 * 1.2) : (8 * 3.5) . (36) 
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It is important to investigate these ratios separately for the production of longitudinally 
polarized vector mesons, where hard physics dominates, and for transversely polarized 
vector mesons, where the interplay between soft and hard physics is most important. 



VII. PRODUCTION OF EXCITED VECTOR MESON STATES 



Equation (0) is applicable also to the production of excited vector meson states if their 
masses Myi satisfy the condition that My, <^ Q^. In this limit, it predicts comparable 
production of excited and ground states. There are no estimates of rjyi] so, following the 
above discussion, we assume that the qq wave functions of a p', uj' and 0' are close to 
the asymptotic value, and, as a rough estimate, we will assume that rjyi = rjy. Using the 
information on the mesons' decay widths from the Review of Particle Properties and 



results of an analysis by Clegg and Donnachie of the properties of excited light vector 
mesons, we find that: 



p(1450) : p° 


^ 0.45 -0.95 , 


(x;(1420) : u 


^ 0.46 , 


p(1700) : p° 


^ 0.22 ±0.05 , 


cu(1600) : u 


^ 0.48 , 


0(1680) : 


^ 0.85 , 


: J/m 


^ 0.5 . 



(37) 



In view of substantial uncertainties in the experimental widths of most of the excited 
states as well as in the values of riy/, these numbers can be considered as good to about a 
factor of two (except for the where Ty is well known). For example, if we would just 
sum over the partial widths of the C(j(1600), as listed in the Review of Particle Properties 
p6| , we would get to'(1600) : u > 0.93. Note, also, that in quark models one would expect 
p(1450) : p ^ u;(1420) : u and p(1700) : p ^ cj(1600) : u, which is not satisfied in the 
cases of the p(1700) and cu(1600). 
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Note that the essential transverse momenta in the qq wave function of the excited 
states are hkely to be significantly larger than in the ground state. This is dictated by 
the spectrality condition in the bound state equation for the minimal Fock qq component 
of the light-cone wave function of the vector meson, which ensures that 

—-f > Ml . (38) 



The estimate which follows from this equation, y (/cf) ^ My//2, shows that the transverse 
momenta characteristic for the minimal Fock component of the wave function should 
increase with the mass of the excited state. Hence, we made a model estimate for the 
dependence of the ratio p(1700) : p taking model 1 type wave functions for the fcj- 



dependence {4>{kt) oc ^^:j_^^2-)2 as introduced in Sect. Ill) with y{k^{p)) = 450 MeV/c, 
and y {ki{p)) / {kf{p')) = Mp/Mpi. The corresponding results for the dependence of 
the ratio p(1700) : p are displayed in Fig. 13 as a solid line. We also present in this figure 
the ratios for \1'' : J/^ and J/^ : p production, normalized to their asymptotic values at 
Q"^ ^ oo given in Eqs. (|36|) and (p7|). 

Larger average momenta in the excited mesons' wave functions imply that the impor- 
tant virtualities in the gluon density essential for a given are larger for the production 
of excited states. This should manifest itself also in an increase of the ratios V : V with 
decreasing x, since the parton distributions depend on {kf{V)), and they increase faster 
at larger (k'flV)). 

For those excited states, which in the nonrelativistic quark model are just radial exci- 
tations of the corresponding ground state (presumably the p(1450) and the ^7(1420)), and 
in the region where one can apply perturbative QCD, our predictions for the increase 



of the V' : V ratios are qualitatively similar to whose of Nemchik et al. ||Tj], who predict 
an increase of p' : p to 1 at large Q^. At the same time, in the calculations of Ref. [|I4|, 
the effects of the quark Fermi motion in the diffractively produced vector mesons, which 
constitute the dominant contribution to the production ratios at larger Q^, were not ac- 
counted for. As a result, the saturation of the V : V ratio occurs in the model of, e.g., 
Ref. already at ~ 10 GeV^. 
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FIG. 13. Relative yields in the diffractive electroproduction of various vector mesons nor- 
malized to their asymptotic values at — > cjo given in Eqs. (|36|) and (l37|) . 

Another important difference stemming from the use of nonrelativistic oscillator wave 
functions, which do not take into account the evolution of the minimal Fock space 
wave functions (which is a strict consequence of QCD [l^), is the qualitative difference 
between the behavior of radial and orbital (e.g. Z^-state) excitations. The production 
of the latter is strongly suppressed at large in the approximation considered in Refs. 
||T3| , p^JT7| , while we find essentially the same rate for both classes of excited states. So 
without mentioning this explictly, the authors of Refs. considered p' and (p' 

corresponding to radial excitations only.0 According to the current interpretation of the 
data, the p(1450) and to'(1420) are radial excitations, while the p(1700) and u;(1600) are 
D-wave states. Note that for the cases where the uncertainties in the analysis of Clegg 

are smaller - u;(1420) and ^^(1600) - we find that the 
are comparable. 



and Donnachie for Ty^e+e 
expected yields at large as given by Eq 



^We are indepted to B.Z. Kopeliovich for clarifying what p'-state they considered. 
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To summarize, a substantial production of excited resonance states is expected at 
large at HERA. A measurement of these reactions may help to better understand the 
dynamics of the diffractive production of vector mesons as well as the light-cone minimal 
Fock state wave functions of these excited states. It may provide the first three dimensional 
images of these states, and it would also allow to look for the second missing excited (f) 
state, which is likely to have a mass of about 1900 MeV to follow the pattern of the p, uj 
and J/\l/ families. 

The relative yield of excited states produced by virtual photons is expected to be 
higher than that for real photons since the Vector Dominance Model (VDM) together 
with Eq. (0) leads to 

airi + N + N) ainl + N + N) 

a{^ + N^V' + N) ai^L + N^V + N) " 

M'y vim - M^y ■ ^^'^ 
In the last step, we used the empirical observation that for the effective cross sections 
of V'N and VN interactions, which enter in the VDM modeLPl < 1, and that 

' 'I — I a-tot(VN) — ' 



^'^Note that these effective cross sections have no direct relation to the genuine interaction 
cross sections. For example, based on geometrical scaling, one expects the interaction cross 
section to increase with the size of the projectile approximately as R = atot{^' / (Ttot{J / ~ 
B^i/B?jj^ ~ 4. However, if one applies the VDM to the extraction of cross sections from 
photoproduction of J /'^ and one finds ~ 0.7 ~ Mj^^/Mg,. This trend seems to reflect 
effects of color screening in the production of heavy quarkonium states |4^Jll. Note, also, that 
photoproduction data do not resolve the p(1430) and p(1700). In the case of p' photoproduction 
off nuclei, similar nuclear absorption effects are observed for the production of the p and p\ 
indicating a{p'N) ~ a{pN). At the same time, application of the VDM to the process 7 p — > p'p 
leads to a{p'N) w (0.37 ± 0.07) a{pN), where we used the analysis of Ref. @ of the V 
e^e~ widths and the data of Chapin et al. |4^ for diffractive p' production. The observed 
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rjy and rjv' are close to their asymptotic values for light mesons while for heavy quark 
systems the values of rjy should be closer to the static limit of r/y = 2. 



Equation (12) is applicable also to vector meson production in weak processes. Con- 
sider, for example, the diffractive production of a D*^ = [cs) meson in W^N scattering. 
To calculate this cross section, it is sufficient to substitute in Eq. (0) the electromagnetic 
coupling constant with g cos 9c, where 9c is the Cabibbo angle. 

It is often discussed that the production of hadrons in e"^e~ annihilation at intermediate 
masses {M^ < 2-^2.5 GeV) can be described as the production and decay of several vector 
meson resonances. Within this approximation, we can generalize Eq. (|T2|) to the case of 
diffractive electroproduction of hadron states in the continuum. To deduce the relevant 
formulae, we assume that for excited vector meson resonances rjy is close to its asymptoic 
value, rjy = 3. This assumption really follows from the application of the dispersion sum 
rule approach. In the limits M| < 5GeV^ Q^, x <^ 1 and oo, we obtain: 

2 

• (40) 



dtdM^ 



t=o AN^aEAiQ^ 
In is worth emphasizing that in the case of diffraction to high enough masses, where 
the intermediate 7* state can be approximated by a free qq pair, the contribution of 
large transverse distances is not suppressed and the cross section gets both soft and hard 
contributions. Consequently, in the intermediate M| range, where resonances are still 
present but the continuum is already essential, one should expect a faster increase with 
energy (at fixed Q^) of the resonance contribution. 



pattern indicates that the production of p' is dominated by average quark-gluon configurations 
(large absorption cross section), while the probability of these transitions is suppressed since the 
transition 7* — > V emphasizes the role of small configurations. Another possible interpretation 
(which maybe complementary to the one above) is a large contribution of nondiagonal transitions 
p' ^ p in the scattering off nuclei. CT arguments (p. 170) suggest that V ^ V transitions 
could indeed be very appreciable. 
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To summarize, the investigation of exclusive diffractive processes appears to be the 
most effective method to measure the minimal Fock state qq component of the wave 
functions of vector mesons as well as the light-cone wave functions of any small mass 
hadron system having angular momentum one. It seems that the cross sections of these 
processes are rather sensitive to the three dimensional distribution of color within the 
wave functions of vector mesons. This observation could be very helpful in expanding 
methods of lattice QCD into the domain of high energy processes. 

VIII. THE KINEMATICAL BOUNDARY FOR THE REGION OF 
APPLICABILITY OF THE LEADING LOGARITHM APPROXIMATIONS TO 

QCD 

In the above, we have explained that the increase of the cross section with energy (^) 
at fixed originates from the energy dependence of the qqT cross section, see Eq. (^) 
and Fig. 3. But, in the region dominated by hard physics, the value of this cross section 
is restricted by the unitarity of the ^-matrix, and even more stringent restrictions follow 
from the condition that the leading twist term should be significantly larger than the next 
to leading twist term. In detail, let us consider the scattering of a small object, a qq 
pair, from a large object, a nucleon. If only hard physics were relevant for the increase 
of the parton distributions at small x, the radius of the gg-nucleon interaction should be 
practically independent on energy, because, in this case, Gribov diffusion in the impact 
parameter space would be determined by the scale characteristic for hard processes. We 
shall explore this observation to establish the region of applicability of the leading In ^ 
and In -r^ — approximations. 

Thus, the geometrical restriction is that, in the hard regime, the inelastic cross section 
cannot exceed the geometrical size of the target nucleon. An equivalent, but practically 
more convenient, formulation of this important property of hard processes is that the slope 
of the differential cross section should not increase with energy. It is just this property 
of hard processes which gives the possibility to deduce an unitarity limit for their cross 
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sections and/or on the region dominated by hard physics. 

We will start our considerations with the derivation of the boundary for the region 
where the decomposition over powers of ^ can be a sensible approximation. A formal 
method to deduce this restriction is to consider scattering of heavy quarkonium QQ states 
from a hadron target, and to vary the mass of the heavy quark. Since the result of such 
a derivation is evident, in this paper we use an equivalent but shorter way to deduce this 
inequahty 0: We establish the limit by applying the optical theorem to calculate the 
elastic cross section for scattering of a small-transverse-size qq pair off a nucleon: 

^^^'(ggiV) = ^e^*(l + /5') , (41) 



dt ' ' 167r 
where 

^ . ,^ . TT (iln(a;(jAr(x, Q^)) , 
13 = ReA,^N/ImAggN ^ ^ ^ rfinV " ^ ^ 

For simplicity, we parametrize the dependence of the differential cross section on the 
momentum transfer t, at small t (which dominate the cross section), in the form of an 
exponential. As the observed t dependence of hard diffractive processes is close to that 
given by the square of the nucleon's electromagnetic form factor, i.e., ~ exp(^t), to 
visualize the deduced limit, it is convenient to express the experimental value of the slope 
of the t dependence of the cross section through the average quadratic radius of a nucleon: 

2 

B ^ 4 ^ 5 GeV^^ ~ ^ where = 0.8 fm. It follows from QCD (cf. discussion in Sect. 
III.A) that this value should be universal for all hard processes. Existing data on hard 
diffractive production of vector mesons are in line with this QCD prediction, and we can 
use it to calculate the elastic cross section: 

crei = ^{l + f3')^^{l + f3'). (43) 
IbirB 47rr^ 

In the derivation of the boundary for the region where the leading In - and In -jS — ap- 
proximations are applicable, one should actually add to Eq. (|^) the contribution from 
proton dissociation, i.e., one should multiply a^i by (1 + 7) where 7 ~ 10 ^ 30%. 
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The necessary condition for the apphcabihty of the decomposition of the deep inelastic 
cross section over powers of is that the leading twist contribution should exceed non- 
leading twist terms. In our case, this implies that Oinei ^ cre«(l + 7)5 where 7 accounts for 
the dissociation of the proton. So, the boundary for the applicability of the decomposition 



we obtain: 



over powers of l/Q^ is = ?/ ^ 1. Using Eq. ( ^31) and above discussed value of -B, 



^^"'^'<<(l+/32)(l + 7)(l + ^)2- ^^^^ 
Thus the region of applicability of the QCD evolution equation is: 

The unitarity of the S'-matrix, which corresponds to the condition ?/ = 1, was discussed 
in Refs. [^,0,^. It can only be translated into an actual kinematical boundary if one 
considers purely hard physics, where B is restricted by the nucleon size. In general, 
however, the slope B and its energy dependence have to be evaluated separately, which 
has not been done so far. Such analysis would lead to the Froissart limit, with B and the 
cross section slowly increasing with energy. Note also that Eq. (|^) has been deduced in the 
approximation where only the leading power of 6^ has been taken into account. A more 
stringent limit follows from the requirement that the single particle density in the center 
of rapidity should be positive. Really, it follows from the application of Abramovskii, 
Gribov and Kancheli (AGK) combinatorics that this multiplicity is proportional to 
c"mei — 4(cre; + CTdiff)- Thus, in the approximation when only one hard rescattering is 
accounted for, we have y < 1/4. 

Therefore, the kinematical boundary for the region of applicability of the leading log- 
arithm approximation is significantly more stringent than that given by the unitarity 
condition: 

Eq. (^61) allows us to estimate down to what x = xum. the equations for ai and ct^^at^pat 
can be applied. For certainty, we choose 7 = 0.2 in our evaluation. Results of a calculation 
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of XiimXQ'^) based on Eq. (^) for the two most current gluon parametrizations are 
presented in Fig. 14. The boundary is even more stringent for p-meson production since 
the average h{Q) are larger in this case (see the dashed and dotted curves in Fig. 14). 
Actually, it is likely that the slow down starts earlier, in particular due to dispersion in h. 
However, a detailed analysis of this effect is beyond the scope of this work. 
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FIG. 14. The boundary of the region of applicability of pQCD for a nucleon target, xnm^Q), 
as determined from Eq. (p6|) for deep inelastic scattering (ul) as well as diffractive electropro- 
duction of p-mesons {"flN — > pN). 



The advantage of the approach used in this paper, as compared to previous attempts 
(for a review and references see Refs. p^ , |55| ), is that we actually deduce the QCD formulae 
for the cross section for scattering of a qq pair off a hadron target. For this cross section, 
the geometrical limit (including numerical coefficients) unambiguously follows from the 
unitarity of the S'-matrix, that is from the geometry of the collision (under the assump- 



tions made in Refs. fS^JS^ on the dominance of hard physics near the unitarity limit). 
As a result, we obtain a more stringent inequality which contains no free parameters. 
Recently, a quantitative estimate of the saturation limit was obtained p6| by considering 
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the GLR model [0,0 of nonlinear effects in the parton evolution and requiring that the 
nonlinear term should be smaller than the linear term. The constraint obtained in that 
paper is numerically much less restrictive than our result: the parameter characterizing 
the scale of corrections found in Ref. [0] differs from ours by a numerical factor ~ jq. 
The qualitative reason for this difference is that the GLR model neglects elastic rescat- 
terings of the qq pair, and it does not take into account the requirement of positiveness 
of the single particle density. The consideration of a specific cut of the double scattering 
diagram, which corresponds to the shadowing of the single hadron multiplicity, and the 
requirement that the single multiplicity should be positive leads, within this model, to 
the restriction that aei + cr^iff cannot exceed l/4crj„e«- The numerical coefficient follows 
then from the application of Abramovskii, Gribov and Kancheli (AGK) combinatorics 
53| . This observation indicates that, in the GLR model hard double scattering 



may decrease the parton density by not more than 25% (and by no more than the ob- 
served diffraction fraction of the cross section which is ~ 15% at HERA). As usual, if the 
corrections become large enough, one cannot restrict the consideration to one particular 
contribution. 

We want to point out that the deduced constraint implies only a restriction on the 
limiting behavior of the cross sections for hard processes, but it gives no clues whether this 
boundary can be achieved. In reality, the increase of the cross sections of hard processes 
should slow down already at larger x than those given by the constraint which we deduced 
here. This is because in the derivation of the boundary we ignore, amongst others, the for 
Feynman diagrams characteristic diffusion in the parton ladder to large distances, effects 
related to the dispersion of b as well as other hard diffractive processes which also fastly 
increase at a; — > 0. The dynamical mechanism responsible for the slowing down of the 
increase of the parton distributions, so that they satisfy, for instance, Eq. (0), is subject 
of discussions. In particular, the triple Pomeron mechanism for shadowing suggested in 



Ref. |]5^ does not lead to large effects at HERA energies, especially if one assumes a 
homogeneous transverse density of gluons ||5^J5^ . 
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IX. LIMITS OF APPLICABILITY OF THE QCD EVOLUTION EQUATIONS 
TO THE SCATTERING OFF HEAVY NUCLEI 



The application of similar consideration to the scattering off nuclei leads to a more 
severe restriction on the gluon distribution in nuclei: 

-'•'"(ts^) £ (1 + 5.) (1 + + ■ (4^) 

Here, Ra is the radius of the nucleus and stems from the dissociation of the nucleus. 
P = ReAgqN / ImAqqN is determined from the x dependence of the gluon density and 
should be universal within a considered approximation. For a heavy nucleus, we may safely 
substitute y = 1 since the total and elastic cross sections are then not far from the black 
body limit, where ~ 2tiR\ and cr^ ~ ^ri?^, even for moderate values of aineiiqqN) . In 
the case of a nuclear target, the limit we consider here is worth a separate investigation 
since interesting effects may arise in the kinematics where the parton densities of a nucleon 
are still far from the unitarity limit. A.Mueller and J.Qiu were the first to draw 
attention to nuclear enhancement of nonlinear effects in the QCD evolution equation. 

The inequality of Eq. (^) corresponds to a rather stringent restriction on the gluon 
distribution in heavy nuclei: 

0^, < g ^ JTT^W^ • '''' 

where Ra = vqA^/^ and tq ~ 1.1 fm. For Q"^ = 10 GeV^ 0:5 = 0.3, 7^ = 0, and y = I, this 
yields 

1 7^ 

jxGa{x, = 10 GeV^) < . (49) 

A more accurate estimate for this constraint could be obtained by using the discussed 
inequality for each b and then averaging over b but not just taking the value in the mean 
point. 

In the above, to quantify the kinematical bounds for and vector meson production, 
we neglected the dispersion over b in the corresponding amplitudes. But, in practice, we 
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need the amplitude in /cj-space, i.e., the Fourier transform of the cross section which we 
discussed above with weights given by the wave functions of the photon or vector meson, 
respectively. For a more accurate estimate, we need to apply the kinematical limit for fixed 
b and then calculate the cross sections based on Eqs. (|T3|) and ( [T^ ) with aggT satisfying 
the bound. 

Let us restrict our analysis here to the case of heavy nuclei as targets, where these 
effects become important at significantly larger x than for a nucleon target, and where 
uncertainties in the estimate of the region of applicability of the evolution equation are 
comparatively small. In fact, we will focus our analysis on the region which would be 
accessible for HERA with nuclear beams and for the LHC: x G 10"'^ ^ 10~^. In this 
kinematics, the gluon distribution in the nucleon is already known reasonably well from 
the recent HERA experiments. 

To estimate nonlinear effects, we allow aqqA{x, h) to reach the unitarity bound, i.e., we 
set 1/ = 1 in Eq. (^), and if cTqg^(x, 6) > vri?^, we substitute it by vri?^. Since in the 
discussed kinematics, the unitarity corrections to CTqqN are still small, the ratio of cross 
sections calculated with and without this unitarity correction characterizes the effective 
gluon shadowing for ctl and vector meson production. In the case of p-meson production. 



we chose the same wave function as in Sect. HI with J (k^) = 450 MeV/c. We also checked 



that varying J (kf) in the range 300 -i- 600 MeV/c practically does not change our results. 



For the sake of an easier comparison of the shadowing effects, in the case of p-meson 



production, we present plots for ^Rv{A,x,Q'^). Clearly this estimate is an upper bound 
for these ratios since it neglects absorbtive effects for the case when the interaction is not 
black. 

To estimate nuclear shadowing within a dynamical model, we also consider the 6-space 
eikonal approximation for the calculation of the same quantities. We first calculate the 
total cross section for the qqN interaction by means of the optical theorem. 





The relevant ratios, x, Q^) = -j^^ and RviA,x,Q'^) = -j2^^^^^^+w^v+ir^ are 

a L ^ da L jdt 

presented in Fig. 15 for A = 40 and A = 238 and for Q"^ =2, 5, 10, 20 and 40 GeV^. 
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tot 



16nB 



{l+P^){l + jA)=(Ttot , 



(50) 



and then evaluate the cross section of the qqA interaction for a given impact parame- 
ter within the 6-space eikonal approximation using the famihar Glauber-Gribov formula 
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Re 



Rl{A,x,Q' 



Rv{A,x,Q' 



A 



A^ 



Jd% dz \ilj^*^(h, z)\'^ (Jtotip) 
Jd% dz {b, z) ijp{b, z) Jd'^B (2 - 2e-(i--/^) '^"''r'^' ^ 



Jd% dz i^^^ib, z) ijjpib, z) atotib)(l - i/S) 



(51) 
(52) 



Here, T(B) = dzp{z, B) is the optical thickness of the nucleus and p{z,B) is the 
nuclear density normalized so that dz d'^B p{z, B) = A. We use Eqs. ( ^1) and ( ^2]) to 
estimate the nuclear shadowing for ai and for p- meson production employing a realistic 
nuclear density distribution, p(r) ~ 1/(1 + exp[(r — i?^)/a]), where Ra = 1.1 A^^^ fm and 
a = 0.56 fm. Corresponding results are presented in Fig. 16. 

One can see from Figs. 15 and 16 that significant gluon shadowing is expected at 



small X, and that shadowing is larger in the p-meson channel, i.e., y i?y(A, x, Q^) < 
Rl{A, X, Q"^). This is again due to the larger b for p-meson production - see the discussion 
in Sect. III. Also, in line with expectations of the color transparency logic, at fixed 
large Q^, shadowing increases with decreasing x - i.e., Rl,Rv <^ 1 for x — and color 
transparency disappears for small enough x - while, for fixed small x, shadowing decreases 
with increasing - i.e., Rl, Rv — > 1 at large and color transparency appears for large 
enough Q^. 

Furthermore, from Figs. 15 and 16, we observe that the nuclear shadowing which we 
deduce from the unitarity limit and from the eikonal model is significantly smaller than 
the black disk limit. This is partly due to the diffuse edge of the nucleus, and a practical 
conclusion from this observation is that nuclear shadowing of the parton distributions 
should be considerably larger for central AA collisions where the effects of the nuclear 
surface should be insignificant. 



56 




X X 
FIG. 16. Nuclear shadowing in DIS in the gluon channel and in diffractive electroproduction 

of vector mesons is evaluated in 6-space eikonal approximation. 
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Note, however, that even those estimates are still likely to underestimate the relevant 
shadowing since they do not explicitly include gluon shadowing in the leading twist, i.e. 
logarithmically decreasing at large . Still, our estimate corresponds to significantly 
larger gluon shadowing than that obtained in a number of previous papers, where it was 
assumed that for ^ few GeV^ the leading twist shadowing in the gluon channel and 
for F2A are approximately the same []62|,p],|63H65[l . Note that qualitative arguments in 
favor of larger shadowing in the gluon channel were presented in Ref . . The quantities 
i?L and Ry of Eqs. (^) and ( ^2|) were also evaluated in Refs. |13,67| using the two- 
gluon-exchange model, in which the elementary cross section of the interactions with the 
constituent quarks does not depend on energy. As a result, in this model, attenuation 
does not depend on x for small x. 

In principle, shadowing effects in the kinematics beyond the region of applicability of 
the evolution equation may change the dependence of the differential cross section on t, 
the momentum transferred to the nucleus In that realm, the eikonal approximation 
leads to the famous pattern with minima and maxima. Another well known feature of 
the eikonal approximation, which is specific for potentials lastly increasing with energy, is 
the shrinking of the diffractive peak with energy. Since, within this model, both soft and 
hard physics are important, the radius of a nucleon may increase with energy, and, for 
this regime, the unitarity restriction for processes dominated by hard physics, which we 
discussed above, could be exceeded and the cross section may increase with x, but more 
slowly as a power of In x. For a nuclear target, such effects are expected to be insignificant, 
although for a nucleon target they may be large. 

To summarize, the analysis presented in this Section demonstrates that large shad- 
owing effects should be present at small x in the gluon channel. This has a number of 
implications for high-energy nucleus-nucleus collisions, where the predictions for hadron 
production strongly depend on gluon shadowing, as was pointed out, for example, in Ref. 
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